Transport criticality at Mott transition in a triangular lattice Hubbard model 
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We study electric transport near the Mott metal-insulator transition in a triangular-lattice Hub- 
bard model at half filling. We calculate optical conductivity (j{u}) based on a cellular dynamical 
mean field theory including vertex corrections inside the cluster. Near the Mott critical end point, 
a Drude analysis in the metallic region suggests that the change in the Drude weight is important 
rather than that in the transport scattering rate for the Mott transition. In the insulating region, 
there emerges an "ingap" peak in o[ui) at low uj near the Mott transition, and this smoothly con- 
nects to the Drude peak in the metallic region with decreasing Coulomb repulsion. We find that the 
weight of these peaks exhibits a power-law behavior upon controlling Coulomb repulsion at the crit- 
ical temperature. The obtained critical exponent suggests that conductivity does not correspond to 
magnetization or energy density of the Ising universality class in contrast to several previous works. 

PACS numbers: 71.27.-|-a, 71.30.-|-h, 72.10.-d 



I. INTRODUCTION 

The Mott metal-insulator transition is one of fasci- 
nating phenomena in strongly correlated electronic sys- 
tems. It is realized in two ways; bandwidth-control 
at fixed electron density and electron-density-control at 
fixed Coulomb repulsion and bare bandwidth. In this pa- 
per, we focus on the former type. The Mott transition 
occurs when the Coulomb repulsion predominates over 
the electronic kinetic energy and it is not accompanied 
by magnetic transition in the original sense. This takes 
place inside the paramagnetic phase at finite temperature 
and magnetic frustration plays an important role for its 
realization. Its phase diagram in the parameter space of 
temperature and (chemical) pressure has insulating and 
metallic phases which are separated by a line of first-order 
Mott transition at low temperature. This phase bound- 
ary terminates at a critical end point, namely, the Mott 
critical end point. This behavior has been observed in 
several materials such as a classic three-dimensional Cr- 
doped V2 0^2 and quasi-two-dimensional K-type organic 
compounds.!^ This is similar to the liquid-gas transition 
in classical liquidJ^ and critical behaviors are expected in 
various properties around the Mott critical end point. 

In the Mott transition, double occupancy (doublon 
density) d = {ni^rii^) plays the role of order parameter, - 
where denotes the electron density operator at site 
i with spin a. Recently, the existence of critical behav- 
iors of double occupancy near the Mott critical end point 
has been reported in some theoretical works. Its critical- 
ity, "Mott criticality," is analyzed with respect to two 
control parameters. One is simply temperature, and the 
other is the conjugate field of the order parameter. In 
our case, the order parameter is double occupancy, and 
this appears in the Hubbard Hamiltonian with multi- 
plied by Coulomb repulsion U . Therefore, U plays the 
role of the conjugate field in the theory of Mott transi- 
tion. In experiments, however, it is not easy to control 
Coulomb interaction directly. It is standard to control 



pressure, instead, since this changes the bandwidth W 
of the material and therefore control effectively the ratio 
U/W. One should note that applying pressure corre- 
sponds to decreasing effective U . A Ginzburg-Landau 
analysis^ or an scalin g an alysis based on the dynamical 
mean field calculations^ have suggested that double oc- 
cupancy shows the same scaling behavior as the Ising 
order parameter and Coulomb repulsion corresponds to 
magnetic field in the Ising model. This is equivalent to 
the liquid-gas transition in classical liquids. However, 
the Mott criticality in other properties is not well un- 
derstood and transport properties are particularly im- 
portant among them. This is the main issue of this pa- 
per and we will report our numerical study on the Mott 
criticality in electronic transport of a triangular-lattice 
Hubbard model. 

To study electric transport in strong correlated elec- 
tronic systems, not only dc-electric conductivity but ac 
response, i.e., optical conductivity, provide useful infor- 
mation on the charge dynamics, in particular, effective 
mass, transport scattering process and electric struc- 
ture. Optical conductivity experiments for several K-type 
organic compounds have examined so far. Insulating 
regular triangular lattice compound k-(ET)2Cu2(CN)3 
shows spin liquid behaviors with no magnetic long-range 
order .I^ISHini Its optical conductivity does not show a clear 
gap but decays smoothly toward zero frequency.!^ In con- 
trast, K-(ET)2Cu[N(CN)2]Cl is less frustrated due to a 
distortion in the triangular structure and exhibits an an- 
tiferromagn etic o rder. Its optical conductivity has a clear 
gap instead.'i2il3l 

Recently, two experimental studies have reported criti- 
cal behaviors in dc-electric conductivity ctq near the Mott 
critical end point, but the two conclusions are not consis- 
tent to each other. Limelette et al. have assumed that uq 
has the same singularity as order parameter in the Mott 
transition and analyzed its singularity in (Vi_xCrx)203. 
Translating the criticality of magnetization with respect 
to temperature and magnetic field in the Ising universal- 
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where t = (T - T*)/T* and p = (P - P*)/P*. 
{aQ,T*,P*) denote the values at the Mott critical end 
point. Recall that pressure plays the role of a con- 
jugate field to (7o in their experiment. The esti- 
mated values of the critical exponents are {S, /3, 7) 
(4.8, 0.33, l-Zp' and they agree with the values for the 
three-dimensional Ising universality class. A similar anal- 
ysis was carried out by Kagawa et al. for the quasi- 
two-dimensional anisotropic triangular lattice compound 
K-(ET)2Cu[N(CN)2]Cl. However, the estimated critical 
exponents are now (5, /3,7) ^ (2, 1, 1),'^^ different from 
the values obtained by Limelette et al. 

Theoretical understanding of these inconsistent results 
is still controversial. Imada et al. argued that critical 
exponents by Kagawa et al. are due to a marginal quan- 
tum critical point in the two-dimensional system .ISEfillZl 
Papanikolaou et al. proposed another theory that the 
critical exponents deviate due to large subleading cor- 
rections related to energy density.l^ The first numerical 
study was performed by Sentef et al. based on the cellu- 
lar dynamical mean field theory and it was claimed that 
the calculated imaginary-time Green's function shows the 
same criticality as the experimental data of 

Among these theories, the first two are semi- 
phenomenological approaches, and in the last one the 
comparison was not made for the identical quantity. 
Therefore, it is important to study the critical exponent 
of conductivity directly by numerical calculations. How- 
ever, it is very difficult to obtain the precise exponent for 
conductivity, since we need to employ some approximate 
scheme like cellular dynamical mean field theory with a 
small cluster size. One important issue is whether the 
transport critical exponent agrees with the one of the or- 
der parameter, and this is the problem that we study in 
this paper. To this end, we calculate these critical ex- 
ponents within the same approximation. Of course, the 
obtained values are not exact critical exponents, due to 
approximation, but if they differ to each other, this im- 
plies that their exact exponents are also distinct. We will 
actually demonstrate that this is the case by using the 
the cellular dynamical mean field approach in this paper. 

To take into account correlation effects and frustrated 
lattice geometry, we calculate optical conductivity based 
on the Kubo formulcP^ and the cellular dynamical mean 
field theory including vertex corrections. This is one of 
the first achievements of numerical approach for conduc- 
tivity, which is a big challenge in numerical computa- 
tions. 

Our paper is organized as follows. In Sec. we first 
review our numerical method, the cellular dynamic mean 
field theory. Then, we explain how to calculate optical 




FIG. 1: (Color online) (a) Isotropic triangular lattice and (b) 
three-site triangular cluster model employing CDMFT. (c) 
First Brillouin zone of an isotropic triangular lattice. Gray 
zone draws the reduced Brillouin zone for the sublattice using 
the three-site triangular cluster. 



conductivity in the cellular dynamical mean filed theory 
including vertex corrections. In Sec. |III[ we discuss ther- 
modynamic criticality near the Mott critical end point. 
In Sec. |IV[ we examine electronic properties near the 
Mott critical end point. In Sec. |Vj we present our re- 
sults of optical conductivity with controlling Coulomb 
repulsion. Sec. |VI| is devoted to the examination of a 
low-energy small peak in optical conductivity on the in- 
sulating side. In Sec. |VIH we present our result of scaling 
analysis of optical conductivity. Sec. |VIII| is summary 
and discussion. 



II. MODEL AND METHOD 

A. Cellular dynamical mean field theory 

The model we consider in this paper is a single-band 
Hubbard Hamiltonian on an isotropic triangular lattice 
at half filling. 
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Here, t is the nearest-neighbor hopping amplitude and 
U is the on-site Coulomb repulsion, /i is the chemical 
potential and is set by tuning electron density to half 
filling in the case of non bipartite system. Cio- is the 
electron annihilation operator at site i with spin cr, and 
nia — cl^Cia-. In what follows, U and T are in units of t. 
The energy dispersion of the kinetic term is given by 
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where k is the wave vector. 

Our primary concern is about Mott transition, and it 
occurs in situations where U is dominant over kinetic en- 
ergy, which is proportional to t. To take into account 
both strong electronic correlations and geometrical frus- 
tration, we use the cellular dynamical mean field the- 
ory (CDMFT)!^ Considering the lattice structure, we 
map the Hamiltonian Q onto a three-site triangular clus- 
ter model coupled to an effective medium Qc,{iLOn) with 
Matsubara frequency ia;„ determined self-consistently as 
shown in Figs, [l] (a) and (b). This self-consistency equa- 
tion is given by 
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where K is the wave vector in the reduced Brillouin zone 
for the sublattice as shown in Fig. [I] (c). t(K) is the 
Fourier transformed hopping matrix for the sublattice 
and Sct(*w„) is the cluster self-energy. We compute the 
cluster Green's function [^^(t)]^^ = ~(7Vcq<j(t)c^^(0)) 
by using the continuous-time quantum Monte Carlo 
(CTQMC) method based on the strong coupling expan- 
sion.lSS Here, a and (3 denote sites in the cluster and r is 
imaginary time. The cluster self-energy is calculated via 
the Dyson equation. 
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Here, Qa{i^n), ^(K), Ga{i^n) and Eo-(«w„) are 3x3 ma- 
trices. In our CDMFT calculations, we consider only 
paramagnetic solution. 

The negative sign problem in the MC calculations is 
one of the important issues, particularly on the frustrated 
lattice. We have succeeded in reducing it by using molec- 
ular orbital bases. Our cluster has C^^ point group sym- 
metry, and electron operators in the molecular orbital 
basis are 
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(Clcr + C2a + Cso-), 
(Cla - C3ct), 
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In these molecular orbital basis, the effective medium 
and the cluster Green's function can be represented by 
diagonal 3x3 matrices. Negative sign problem is not 
serious when calculations are performed with these basis 
functions. Figures [2] (a) and (b) show the average sign. 



{sgn) = 
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in the MC calculations on U-T parameter space, where 
Pi is the sign of the weight for MC configurations. One 
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FIG. 2: (Color online) (a) T-dependence of the average sign 
(sgn) for two U's and (b) fZ-dependence of (sgn) at T — 0.10. 
(c) Local Green's function GI\t) at T = 0.06 for U = 7.5 
and 11. Error bars are shown at two typical values of r: r = 
and 13/2. 



can see that (sgn) is large even for smaller U and lower 
T, and this guarantees higher-accuracy calculations for 
various quantities. For example, we have calculated elec- 
tron density n by performing 2 x 10'' MC sweeps and 
averaging over 64 samples. The error An is small and 
the typical value is An/n ~ 0.003 for U = 7.5 and 11 
at T = 0.06. In Fig. [2] (c), we show the diagonal part 
of the cluster Green's function G';^^(t) as a function of 
imaginary time r, for two values of U. Its error at t = 
is {AGI^/GI^){t = 0) - 0.02, and this is larger than the 
one in the direct measurement An/n ~ 0.003 but is still 
quite small. To determine the chemical potential /i, we 
have used the directly calculated value of n rather than 
calculating from G^^. 

After the self-consistency loop of CDMFT converges. 
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FIG. 3: (Color online) Particle current in an isotropic trian- 
gular lattice. 



we calculate the lattice Green's functions Go-(k, iw„) of 
the triangular lattice from the cluster quantities. There 
are several methods for this interpolation, and we use the 
cumulant method in this study. It has been shown that 
the cumulant interpolation works well in a wide range of 
U from metallic to insulating state.'^SHlI! We first intro- 
duce and calculate the cluster cumulant 



(10) 



and proceed to calculate the cumulant on the original 
lattice system 

M,(k,ic.„) = ^^il4"'5(ic.„)e-''"('-°-'-^), (11) 
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where Nc{— 3) is the cluster size and is the real-space 
position of site a. Then, the lattice Green's function is 
given as 



Ga(k,iw„) 
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B. Optical conductivity including vertex 
corrections 

Let us discuss vertex corrections in details. In the 
single-site DMFT, the current-current correlation func- 
tion has no vertex correction and conductivity is calcu- 
lated by simple conv olution of the lattice Green's func- 
tion Go-(k, «w„).l^^'^This is because current is a site off- 
diagonal operator, and this point is different from the 
case of, for example, spin susceptibility, where spin den- 
sity is site diagonal. Effects of vertex corrections on spin 
susceptibility have been studied already in the single- 
site DMFT approach!^-**^^-' Vertex corrections have finite 
effects on co nduct ivity in the CDMFT except in some 
special cases ,'^^'21] b^t in some preceding works of the 
CDMFT, cond uctivi ty has been calculated without ver- 
tex corrections.'^^'^ll Their effects were studied particu- 
larly for the square lattice cases using dynamical clus- 
ter approximation (DCA), which is a cluster extension 



of DMFT in the k-space.l^^El To take into account the 
correlation effects in current response further, we use 
CDMFT and calculate optical conductivity including ver- 
tex corrections. This is a big challenge in numerical com- 
putations. 

We start with explaining our algorithm including the 
vertex corrections in CDMFT. Optical conductivity ten- 
sor cryj/'(a;) is generally defined by the Kubo formulaP^l 
as 
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where e is the elementary charge, uj is the frequency, and 
q is wave vector. Here, Xi^i/' (Qi ^) is the current-current 
correlation function, and is calculated from Fourier trans- 
form of the correlation function in real time t 
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where Jiy(q, t) is the particle current operator in the 
X, y) direction and is defined as follows. 
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where a is the lattice unit and N is the total number 
of sites. r=(a,0), (a/2,V3a/2), or (a/2, -^/3a/2). In a 
triangular lattice system, the particle current is defined 
on each of the three bonds as shown in Fig. [3] and the x- 
component of the net current is given as Ji -I- ( J2 + Js)/^. 
In our calculations, we only focus on a^xii-^) because of 
CTj,^/ (w) = a{uj)5^,y> in an isotropic triangular lattice. We 
denote (t{ui) = axx^i-^) and x(*f^n) = Xxxii^n) later. From 
Eqs. ^^ and (fl^, xit) = x(q ^ 0,t) is given by 
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where is the x-component of velocity described as 
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In order to calculate cr(w), we define the current- 
current correlation functions in Matsubara space a^^ 
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where 
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FIG. 4: (Color online) Diagrams for (a) the current-current 
correlation function x(ivn) in Eq. (181 and (b) the full vertex 
function V^^.iiv^) in Eq. Q. 



In our calculation as shown in Fig.|4](a), Eq.(18) is exact 
except that the dependence on internal frequencies in the 
full vertex function F^^, (w„; iuii,iuj'^) has been averaged 
over, namely, T^^,{iiy„;iuji,iuj'^) w r^^'(w„). The pre- 
ceding studies on the vertex corrections have reported 
that in the case of spin susceptibility, this approxima- 
tion qualitatively reproduced the results consistent with 
full calculation .I^SMl] Since full calculation for the vertex 
function also requires much longer computational time, 
we use this approximation in this study. This is an uncon- 
trolled approximation and it is very important to check 
its validity for the current vertex and this should be stud- 
ied in future. 

To obtain T^^, (w„), we first calculate 
two-electron Green's functions in the cluster 
xTfi^s{r) = {cL{r)cpAr)cl^,{0)c5a'm with imagi- 
nary time t by using CTQMC method, where a-6 
denote sites in the cluster. We then evaluate the 
irreducible vertex function in the cluster lap-ysi'^^n) via 
the Bethe-Salpeter equation, 

aa' I ■ \ O.cra' / ■ \ , \ ^ O.ctct" / - \ 

a' ,7'<5' 

X i^:'f:;s,{z,,„)x;:szis{i'^,,). (20) 
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FIG. 5: (Color online) (a) Current-current correlation func- 
tion xi^^n) and (b) optical conductivity ^{ijj) for two U's at 
T = 0.10. 



then calculate the lattice irreducible vertex 



and obtain r?^,(w„) via the Bethe-Salpeter equation as 
shown in Fig.]4](b), 



k",cr''' ii-Ji 



(22) 



Once r^^, (w„) is obtained, we calculate a{ivn) from 
Eq. ( |l8l ). Finally, we calculate real frequency quantity 
(j{ui) via analytic continuation w„ -> cj -I- zO by using the 
maximum entropy method.'^ In what follows, we nor- 
malize alto) by the unit of /h, where h is the reduced 
Planck's constant, and lo are in units of t. 

In Fig. [5] (a), we show the results of x(wn) for two 
L/'s at r = 0.10. As mentioned in the previous subsec- 
tion, negative sign problem in our MC calculations is not 
serious, and then we can calculate x(if„) with high ac- 
curacy. For example, at ?7 = 8.5, we perform ~ 10^ MC 
sweeps and averaging over 512 samples and the relative 
error Ax(«t'n) is small, e.g. Ax(0)/x(0) '-^ 0.03. By us- 
ing the maximum entropy method, we obtain a{ui) from 
x{ivn) as shown in Fig. [5](b). In Sec.|vj we will explain 
the results in detail. 
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III. THERMODYNAMIC CRITICALITY 

In this section, we will examine Mott criticality in 
thermodynamics by our CDMFT calculations before dis- 
cussing electric transport. 

We first start with identifying the location of the finite- 
T Mott transition and determine U-T phase diagramPH 
To this end, we investigate the singularity of the double 
occupancy 

d = (rii^nii) (23) 

with varying U at several T's. Figure. |6] (a) presents U- 
dependence of d for various T's. At lower temperature 
T = 0.09, d shows a jump and hysteresis, corresponding 
to the first-order Mott transition. At higher temperature 
T = 0.12, d changes smoothly and there is no hysteresis 
in the {/-dependence. This indicates a crossover between 
metal and insulator. We have checked that the jump and 
hysteresis of d shrink with increasing T and disappear 
for T > 0.10. Repeating calculations for various T's, the 
Mott critical end point is found to be at U* = 9.4 and 
T* = 0.10. We determine the U-T phase diagram and 
the result is shown in Fig. [6] (b) . In the following, we will 
analyze a singularity of d around the Mott critical end 
point. 

In the Mott transition, d plays the role of order 
parameter.El Previous single-site DMFT and four-site 
CDMFT calculations show that [/-dependence of d is 
well fitted by a power-law scaling function with the mean 
filed exponent 1/i5mf = 1/3 near the Mott critical end 
point This confirms that the Mott transition belongs 
to Ising universality class. These approaches are still 
a mean field approximation, and therefore the obtained 
critical exponent is the mean field value. Let us examine 
our results of d at T = 0.10. d changes continuously with 
U but exhibits a divergent slope near U — 9.4. We try 
to fit the curve with a power-law scaling function. 



\d-d*\ A±|C/-[/*| 



(24) 



where the subscript — (+) denotes the metallic (insulat- 
ing) side. Figure [6] (c) shows the result of scaling analysis 
for d using five points on each side around the critical 
point. This shows that the fitting is successful on both 
sides, namely, d exhibits the critical behavior. The ob- 
tained critical exponents are 

(l/<5", 1/^+) - (0.32 ± 0.05, 0.30 ± 0.04), (25) 

and our CDMFT calculations also reproduce the mean 
field exponent within accuracy. Here, we estimated 
and those standard deviation by repeating scaling anal- 
ysis for each of the 64 data sets of d. 



IV. ELECTRONIC PROPERTIES NEAR THE 
MOTT CRITICAL END POINT 

In the previous section, we have examined the finite- 
T Mott transition and obtained the U-T phase diagram 
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FIG. 6: (Color online) (a) Dependence on U of double occu- 
pancy d for various T's. (b) U-T phase diagram. Diamonds 
represent the boundaries of the region where metallic and in- 
sulating solutions coexist. Line of the first order transition is 
a guide for eyes. Circles show the crossover above the Mott 
critical end point, (c) Scaling analysis of d upon controlling 
U at T — 0.10. Symbols show calculated data and lines are 
the results of the fitting Eq. p4| . 



as shown in Fig. |6] (b). Next, we will study electronic 
properties, particularly near the Mott critical end point. 
This is useful for understanding electric transport, which 
will be discussed in Sec. |Vl 

Density of states is the most fundamental quantity of 
electronic structure and it has been calculated by sev- 
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FIG. 7: (Color online) Density of states A{u)) for various C/'s 
at T = 0.10. 

eral groups for the triangular lattice.'^^HIlIl These results 
demonstrated the existence of a quasiparticle peak near 
the Fermi level in the metallic region and its disappear- 
ance at the metal-insulator transition resulting into the 
formation of a gap at the Fermi level. However, detailed 
behaviors near the Mott critical end point are not well 
understood, and we examine them. 

Density of states A{ljj) is calculated from the lattice 
Green's function in Eq. as 

^H = -^^ImG(k,c^ + iO). (26) 

k 

Here, the Green's function Go-(k,w) does not depend 
on spin a since we treat only paramagnetic solution in 
CDMFT calculations and we drop the spin index. To 
obtain A{ijj), we have carried out analytic continuation 
iujn uj + iOhy using the maximum entropy method.!^ 

Figure [t] presents the change of A{ijj) with varying U at 
T = 0.10 fixed. At all the values of U , A{(jS) shows broad 
peaks at [/ ^ ±4, and they correspond to the upper and 
lower Hubbard bands. At [/ = 9.3, in addition to them, 
A(uj) shows a sharp peak ai uj — 0. This is a quasiparticle 
peak: indication of the metallic state. At larger U — 10, 
the quasiparticle peak disappears and there emerges a 
gap around u — Q and its existence is a characteristics of 
the insulating state. 

The interesting characteristic is that two small peaks 
appear around w = on the insulating side near the 
Mott transition, for example, dX U — 9.41. With ap- 
proaching the Mott transition, their peaks shift towards 
w = and grow their intensity. These peaks continue to 
the quasiparticle peak on the metallic side. We have also 
confirmed this behavior for various T's, and have found 
that the emergence of two small peaks is clearer near the 
Mott critical end point. This is more complicated than 
the canonical picture of Mott transition. The canoni- 
cal picture is that a quasiparticle peak disappears at the 
metal-to-insulator transition and a gap opens continu- 
ously. Our calculations show that this picture is supple- 
mented by the appearance of the low-energy small peaks 
on the insulating side near the Mott transition. 




FIG. 8: (Color online) Optical conductivity for various 

f/'s at (a) T = 0.12, (b) T = 0.10, and (c) T = 0.09. 



V. TRANSPORT PROPERTIES 

In this section, we will investigate optical conductivity 
a{uj). 

First, let us summarize the general trends of cr(w) in 
the U-T space. Figures [s] (a)-(c) show the variation of 
<j{uj) upon controlling U at several T's. First, metallic 
state is realized for small U at all T's. cr{uj) shows a 
Drude peak at low w as expected, indicating the forma- 
tion of coherent quasiparticle. There exists a broad peak 
aX w ^ U and this comes from excitations to the Hubbard 
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FIG. 9: (Color online) Dependence on U of a Drude weight 
Do (upper panel) and a transport scattering rate 1/r (lower 
panel) at various T's. 



band. Next, insulating state is realized for large U at all 
T's. The Drude peak disappears and its absence is a char- 
acteristic of insulator. The most important part is near 
the Mott critical end point T ~ 0.10 and U ^ 9.4 and we 
investigate <t{uj) there upon controlling U at various T's. 
At T = 0.10, we find that a{uj) at low lu changes quickly 
at ?7 ~ 9.4, At the lower temperature T = 0.09, cr(w) 
at low uj shows a jump, corresponding to the first-order 
transition, and we also confirm that there is a hystere- 
sis. At the higher temperature T = 0.12, 17(0;) at low 
UJ shows a smooth crossover from metallic to insulating 
state. In the following, we will analyze criuj) on metallic 
and insulating sides in detail, particularly, focusing on its 
dependence on f7 at T = 0.10 near the Mott critical end 
point. 



Metallic side 




(oj-co , )/W~ 

^ peak' 



(b) 
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FIG. 10: (Color online) (a) Fitting results of an incherent 
broad peak on the insulating side by Eq. (28 1, (b) Ingap 



peak (JiG(t<j) on the insulating side. At U = 9.4, a small 
Drude peak appears, accompanied by an ingap peak around 
w ~ 1. (c) Dependence on U of the weight of an ingap peak 
Dig for various T's. 



For simple metals at low temperature, it is well known 
that the low-cj part of (j{uj) is described by a simple 
Drude formula 



<j{uj) = Re 



Da 



(27) 



where Dq is a Drude weight and 1/r is the transport 
scattering rate. We fit the low-w peak of our data by 
Eq. ( p7| and examine the changes in Dq and 1/t upon 
controlling U a,t T fixed. 

Figures |9] (a) and (b) show {/-dependence of Dq and 
1/r for various T's. One may expect that with approach- 



ing the Mott critical end point, 1/r increases with in- 
creasing U and diverge. One can see that 1/r does not 
change noticeably with U and shows small increase near 
the Mott transition point. Dq shows noticeable decrease 
with U and drops drastically near the Mott transition 
point at T = 0.10. These results suggest that Mott tran- 
sition is driven by the change in Dq rather than 1/r. 



B. Insulating side 

Now, let us study the insulating side near the Mott 
transition. The most important difference from conven- 
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FIG. 11: (Color online) k-resolved spectral function Ak(i^) 
at three values of (7 at T = 0.10 along the path T-K-M-r in 
Fig. 0(c). 



tional insulators is that (7{lj) has two parts; one is a broad 
Hubbard peak around co ^ U and the other is a small 
peak at low uj. The former peak persists from the metal- 
lic side and comes from the contribution of the Hubbard 
band, (Thb(^^)- Its peak position shifts from the strong 
coupling limit a; ~ C/ by a finite amount ~ —3, which 
conies from the kinetic energy of a pair of doublon and 
holon. We find that (Thb(i^) is well fitted by a simple 
Gaussian form, 



CTHB (^) 
"■(l^poak) 



— Aexp 



-B{w - Wpeak)' 



(28) 



particularly for lo below the peak position Wpeai? of 
fHEl"^)- This is shown in Fig. [lO] (a). Here, A and B 
are constant and W- is the half width below Wpoak- 

In order to carry out detailed analysis of the small 
peak at low w, we subtract (Thb(^) from the total cr(w) 
by Eq. (281 and the remaining part is defined as aiQ(uj). 



peak below w ~ 2. The peak position shifts toward lower 
UJ with approaching the Mott transition point. We will 
call this structure an ingap peak hereafter. For this part, 
we define its intensity by. 



D 



IG 



dujaiQ (lu) . 



(29) 



We analyze our data in the insulating phase and show 
the t/-dependence of Dig in Fig- 10 (c). Dig increases 



with decreasing U and its growth is drastic at T = 0.10 
near the Mott critical end point. 



VI. INGAP PEAK AND SPECTRAL FUNCTION 

We have shown that optical conductivity <j(uj) exhibits 
an ingap peak, an interesting structure on the insulating 
side near the Mott transition. In this section, we briefly 
discuss the ingap peak from a viewpoint of electronic 
structure. 



In Sec. IV we examined the density of states A{uj) 



and found that A{uj) exhibits low-energy two small peaks 
upon approaching the Mott transition point from the in- 
sulating side. To see its origin, let us calculate the k- 
resolved spectral function Ak{uj). Two previous stud- 
ies also calculated ^k('^) in a- triangular lattice and re- 
ported two features with compared to the square lat- 
tice case. One is weaker intensity of quasiparticle peaks 
in the metallic phase, and the o ther is the absence of 
anomolous pseudogap phase .I^SBSl Ohashi et al. studied 
an anisotropic triangular lattice system and observed the 
split of the quasiparticle peak in the insulating phase, 
which is attributed to enhanced magnetic fluctuations.'^ 

Figure 11 presents Ak{uj) at three values of U at 
T = 0.10. The data are plotted along the path F-K- 
M-F in the Brillouin zone shown in Fig. [l] (c). At the 
weakest correlation U = 9.4, in addition to broad peaks 
corresponding to the upper and lower Hubbard bands, a 
quasiparticle peak exists around uj — 0. The energy dis- 
persion is strongly renormalized to a smaller width due 
to strong correlation efl[ects. We have confirmed that the 
Fermi surface is nearly spherical. The quasiparticle peak 
disappears at the stronger correlation U = 9.41 and there 
emerges two small peaks centered around uj — 0, which 
exist in the whole Brillouin zone. At the strongest cor- 
relation U = 10, these small peaks disappear and there 
emerges a clear gap around cj = 0. 

We have found that the small peaks start to appear at 
the same time of emergence of the ingap peak in cr(w). 
We have also checked that the energy difference between 
the two small peaks agrees with the position of the ingap 
peak. Therefore, the ingap peak is an optical transition 
between these two peaks, which is associated with the 



The results are shown in Fig. 10 



(b). They all show a collapse of quasiparticles. 
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FIG. 12: (Color online) Dependence on U of (a) dc-electric 
conductivity (tq and optical weight D for various T's. 



VII. TRANSPORT CRITICALITY 

In this section, we provide a detailed analysis of optical 
conductivity cr(w) near the Mott critical end point to 
discuss the Mott criticality in electric transport. 

We first examine whether a{u)) exhibits a singularity 
at the Mott critical end point. The most important issue 
is how to identify a singularity of cr(aj), and this has not 
been well established. Since <j{ijj) is not derived from free 
energy, protocol of scaling analysis is not uniquely deter- 
mined like that for thermodynamic quantities. Therefore, 
we choose several quantities in (t{lo) as a scaling variable 
and examine their criticality. 

The simplest choice as a scaling variable is dc-electric 
conductivity ctq = = 0). This has been used in 
experime ntal s tudies on the Mott criticality in electric 
transport.'^^'^ The experimental results have demon- 
strated the existence of a singularity of ctq near the Mott 
critical end point. However, uq is not necessarily the best 
choice because impurity scatterings may contaminate its 
intrinsic critical behavior. In fact, scaling analysis of the 
experimental data has been performed only on the metal- 
lic side. Here, our results of J7-dependence of do at vari- 
ous T's are plotted in Fig. Il2] (a). At T = 0.10, ctq shows 
a singularity near U ~ 9.4. However, its scaling analysis 
on the insulating side leads to rather large relative errors, 
since ctq is very small, as we will explain later. 
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FIG. 13: (Color online) Scaling analysis of optical weight D. 
Symbols sh ow calculated data and lines are the results of the 
fitting Eq. (|30l. 



An alternative scaling variable may be the optical 
weight of a Drude peak on the metallic side. This 
quantity is obtained by integrating the low-w peak of 
(t(w), and then is more robust against impurity scatter- 
ings. However, the insulating side does not have the 
Drude peak. In Sec. |V] we have shown that with in- 
creasing t/, the Drude peak on the metallic side changes 
to the ingap peak on the insulating side. Therefore, a 
candidate on the insulating side is the optical weight of 
an ingap peak D\q. We show the ^/-dependence of the 
total optical weight D = Dp + T'ig for various T's in 
Fig. 12 (b). At higher temperature T = 0.12, D shows a 



smooth crossover from the metallic to the insulating side. 
At lower temperature T = 0.09, D shows a jump near 
U — 9.4. We have also checked the existence of hystere- 
sis, and these results indicate the first-order Mott tran- 
sition. At T = 0.10, D changes drastically near V = 9.4 
without showing its hysteresis similar to o-q. 

In the following, we will examine a singularity of D in 
detail, particularly its dependence on f/ at T = 0.10. 



A. Preliminary scaling analysis and accuracy of 
optical conductivity 

Let us start a scaling analysis of the optical weight at 
T — 0.10. To perform more elaborate scaling analysis, 
we estimated error in optical conductivity. In this sub- 
section, we show these results and proceed to the scaling 
analysis of its optical conductivity in the next subsection. 

As mentioned in Sec. |III[ criticality in thermodynamics 
has been investigated by examining the thermal average 
double occupancy c? as a function of T or U . Following 
this analysis, we try to fit the curve of the optical weight 
D with a power-law scaling function, 

\D -D*\= A^\l] -U*\^l^^ , (30) 
where the subscript — (+) denotes the metallic (insu- 
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FIG. 14: (Color online) Optical conductivity for two U's at 
T — 0.10; a{u}) (♦) and d-{uj) {()). Error bars are estimated 
by the jackknife analysis. 



lating) side. D* denotes the value at the Mott crit- 
ical end point {T*,U*) and is a constant. The 
critical exponents l/S^ reflect the underlying universal- 
ity class. Fitting parameters are l/Sf, A±, D* and 
U* . Figure [13] shows this result of scaling analysis for 
Eq.(30l using five points on each side of the Mott crit- 



ical end point. The fitting is successful on both sides, 
and D exhibits a critical power-law behavior. The crit- 
ical region is about \U — U*\ < 0.2, which is similar to 



that of d discussed in Sec. Ill In this scaling analysis, 
the critical end point is determined simultaneously and 
{U*,D*) - (9.40,0.11). The obtained critical exponents 
are {1/S-,1/S+) - (0.15,0.13). 

This is a preliminary scaling and we need error analy- 
sis of D for more reliable value of the critical exponent. 
In the maximum entropy method, the error in imaginary- 
time QMC samplings does not transfer to the calculating 
real-frequency data. Here, we use the jackknife method 
for an error analysis of optical conductivity.'^ In the jack- 
knife method, we divide total N sets of data of a quantity 
O into Nb bins of the bin size m = N/Ni,. For each bin 6, 
we first define the average Ob by taking away bin b from 
the whole data. For these Ni, data, we then define their 
average O and error AO, 



O = 



AO = 



— To 



b=l 



b=l 



(31) 
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FIG. 15: (Color online) Bin size m-dependence of jackknife 
error in optical weight AD. 



nearly constant, one can use that value as the error of O. 
For optical conductivity, we have performed this analysis 
with varying m from 2 to 64 for N = 512 data sets. We 
have checked the dependence of AO on m, and found 
that the errors are almost constant for all w's over the 
whole range of bin sizes. In the rest of the paper, we will 
show the average and the error of optical conductivity 
and other quantities that are analyzed with the bin size 
m = 16 unless explicitly mentioned. 

Figure [14] presents the results of optical conductivity 
CT(a;) at T = 0.10 on the metallic {U = 9.35) and insulat- 
ing [U — 9.45) sides. We find that fT(aj) agrees with cr{uj) 
within accuracy in both cases. The relative error is small 
at high w's, e.g. Act/ct - 0.01 at w ~ C/ = 9.35 and 9.45, 
while the error at small w's is larger, e.g. Aaja ^ 0.12 
and 0.22 at a; ~ for [/ = 9.35 and 9.45, respectively. 
Having the error in optical conductivity, we proceed to 
estimating the error in optical weight and its scaling anal- 
ysis. 



B. Critical exponent of optical weight 

In the previous subsection, we have estimated the error 
in optical conductivity. We will evaluate the error in the 
critical exponents of optical weight, and proceed to their 
discussion in detail. 

From the N = 512 sets of optical conductivity, we 
calculate the optical weight for all and then determine 
their average Z) and error AD by the jackknife method. 
The bin size dependence of AD is shown in Fig. 15 for 
two C/'s. This dependence is nearly constant for this 
and all the other values of C/, and this shows that the 
jackknife error analysis works well in our calculations. 
It is important that the relative error AD / D is much 
smaller than Aa/a, e.g. AD/ D - 0.03 for U = 9.35 
and 9.42. The aver age D agrees with D calculated at 
the beginning of Sec VII within the error AD in Fig. 16 
(a). Note that D is an average defined differently. For 
this, the average is taken over the imaginary-time data, 
and this averaged result is analytically continued to real- 
frequency and finally the optical weight D is calculated. 
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FIG. 17: (Color online) Scaling analysis of (a) dc-electric con- 
ductivity (To and (b) Drude weight Dq. Symbols show calcu- 
lated data and lines are the fitting results. 
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FIG. 16: (Color online) (a) [/-dependence of optical weight; 
D (♦) and D (O)- (b) Scaling analysis of D when errors are 
evaluated. Open symbols show calculated data, dotted lines 
are the fitting results and error bars are estimated by the 
jackknife analysis, (c) Error in the critical exponents 1/S± 
due to fitting. 



Now, with these data of D, let us perform the scaling 
analysis. For each of the Nt, data of the optical weight, 
we repeat the scaUng analysis Eq.pO|), and estimate the 
critical exponents by the jackknife method. Their 



average and error are 

il/S-,l/St) - (0.15 ± 0.04, 0.16 ± 0.03). 

Fig. [l6] (b) shows the fitting with these values, and this 
shows that the fitting is successful on both of the metallic 
and insulating sides. 

Finally, we discuss the value of 



III 



we 



In Sec 

have examined the critical behavior of double occupancy 
and have confirmed the mean-field value of Ising critical 
exponent 1/6 = 1/3. If conductivity has the same criti- 
cal behaviors as claimed by Limelette et al.^ the optical 
weight should have the same critical exponent and there- 
fore 1/Sf = 1/3 in our case. This is because these two 
are calculated in the same approximation, i.e. CDMFT. 
However, our result of the critical exponent is not close 
to 1/3 or any of these values in the Ising universality 
class (1/15 and 1/4.8 for the two- and three-dimensional 
casePl) . 

We have also estimated the error in the critical expo- 
nents due to fitting. Figure 16 (c) shows the accuracy of 
our fitting. The residue of fitting is plotted in the space 
of the two fitting parameters The order parame- 

ters are all optimized for each pair of 1 /Sf . The residue 
is shown by its relative difference compared to the min- 
imum value: i.e. 5% means that the residue is larger by 
5% than that of the globally optimized set of the fitting 
parameters. This plot shows that the exponent of double 
occupancy 1/6 = 1/3 has a very large residue and this 
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value is very unlikely. Thus we can exclude the possibil- 
ity that the optical weight has the same critical behavior 
as the order parameter in the Mott transition. 

Recently, Papanikolaou et al. suggested a different sce- 
nario for the scaling law of conductivity; the subleading 
term due to energy density dominates rather than the 
leading term (magnetization) .1^ If that is the case, one 
may see that the mean-field critical exponent of energy 
density 1/(Se — 1/1.5. However, the obtained critical ex- 
ponent is not close to 1/Se either. 

We have also tried the same analysis of other quantities 
of a{Lu). Figure 17 (a) shows the result for dc-electric 
conductivity ctq. We confirm that ctq shows a critical 
power-law behavior on each of the metallic and insulating 
sides. However, the obtained two critical exponents are 
clearly distinct, {1/Sj - 0.16) < (1/(5+ - 0.07), while 
1/5^ is close to the one in the optical weight. We also 
analyze a Drude weight Dq on the metallic side only by 
the power-law scaling function of U, 



Do = A-{U* - Uf'^^. 



(32) 



As shown in Fig. 17 (b), Dq shows a critical power-law 
behavior and leads to a similar value 1/5'^ ^ 0.13 to that 
in the optical weight. 



VIII. SUMMARY AND DISCUSSION 

In this paper, we have studied the Mott criticality in 
electric transport of the half-filled Hubbard model on a 
triangular lattice. We have succeeded in the numerical 
calculations of optical conductivity taking into account 
correlation effects and frustrated lattice geometry by the 
cellular dynamic mean field theory including vertex cor- 
rections. 

Before calculating conductivity, we checked thermody- 
namic properties and determined the phase diagram and 
the Mott critical end point within our version of numer- 
ical method. Wc then examined the criticality of double 
occupancy, which is the order parameter in Mott tran- 
sition. To discuss transport properties, we also inves- 
tigated electronic structure. In addition to density of 
states, which is a local quantity, we calculated single- 
electron spectra at each wave vector. 

The phase diagram determined in our work is consis- 
tent with the previous works, and the position of the 
critical end point is confirmed. Our calculation also re- 
produced the expected critical exponent of double occu- 
pancy, and that is the mean-field value of the Ising order 
parameter with respect to its dual field 1/(5mf = 1/3. 

The density of states shows two small peaks around 
zero frequency on the insulating side near the Mott tran- 
sition, which originate from a quasiparticle peak on the 
metallic side. This suggests that the canonical picture 
of Mott transition is supplemented by the appearance of 
these small peaks. 

We have extensively analyzed optical conductivity and 
the criticality of its characteristic quantities. In this pa- 



per, we mainly focused on the criticality upon controlling 
Coulomb interaction at the critical temperature. Before 
performing detailed analysis of transport criticality, we 
have tried a Drude analysis of optical conductivity. It in- 
dicates that Mott transition is primarily determined by 
the change in the Drude weight rather than transport 
scattering rate. Another important point is that optical 
conductivity on the insulating side has a small peak at 
low frequencies, which we named an "ingap" peak. The 
wave- vector resolved single-electron spectra show the in- 
gap peak, which originates in transitions between two 
low-energy small peaks around zero frequency appearing 
after the collapse of quasiparticles. 

The Drude peak on the metallic side continues to the 
ingap peak on the insulating side with Coulomb repul- 
sion. We have analyzed the singularity of the optical 
weight of these peaks near the Mott critical end point and 
have compared its critical exponent with that of double 
occupancy within the same calculations. However, the 
obtained critical exponent of the optical weight differs 
from that of double occupancy or corresponding value of 
Ising exponent in any dimension. 

We have justified the belief that double occupancy has 
the same scaling behavior as the Ising order parame- 
ter. The most important point is that the critical ex- 
ponent in conductivity is unconventional in contract to 
that of double occupancy within our numerical accuracy. 
This discrepancy may disprove the appealing and pop- 
ular working hypothesis that conductivity exhibits the 
same scaling behavior as order parameter in the Mott 
transition. Regarding another possibility that conductiv- 
ity corresponds to energy density of the Ising universal- 
ity class, the corresponding mean field exponent differs 
from our result too. 

In the cellular dynamical mean field calculations, the 
cluster size is one of the important conditions. As for 
double occupancy, preceding studies with single- and 
four-site clusters both concluded the same critical ex- 
ponent as our calculations with a three-site cluster.!^ 
Therefore, one may expect that the critical exponent in 
conductivity is also robust against the change of the clus- 
ter size, but we will need to check this point directly using 
different cluster sizes and geometries. 

Another point is about model Hamiltonian. Our re- 
sults are obtained for the Hubbard model, i.e., the stan- 
dard model of a strongly correlated electronic system. 
However, some of the terms neglected in the Hubbard 
model may become important and crucial near a metal- 
insulator transition. They include electron-phonon inter- 
actions and long-range part of Coulomb interaction. If 
they dominate criticality of the transition, the observed 
exponents differ from those of the Hubbard model. Our 
understanding of these models is very limited, and this 
is a future problem. 

Lastly, we make a comment on scaling analysis in ex- 
perimental works. Our scaling analysis of dc-electric con- 
ductivity shows that two critical exponents on the metal- 
lic and insulating sides do not agree. In the experimental 
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studiesp'i^ a singularity of dc-electric conductivity has 
been analyzed only on the metallic side. It is interesting 
to check whether both critical exponents are identical in 
experiments. In optical conductivity, it is expected that 
a low-energy peak appears on the insulating side like an 
ingap peak. It is also interesting to observe this peak and 
examine if it evolves into the Drude peak with varying 
pressure. 
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